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A commercially viable GaAs device technology for field- effect tran- 
sistors, integrated circuits, and lasers is critically dependent on the 
availability of high-quality, single-crystal boules with controlled di- 
ameter. We have modeled a diameter control scheme based on the 
monitoring of crystal weight for liquid-encapsulated Czochralski 
(LEC) growth of GaAs. The presence of the B 2 3 liquid encapsulant 
and significant capillary forces make the direct interpretation of the 
weight-gain signal and its time derivative (DWGS) more complicated 
in comparison with pulled materials such as oxide crystals. We have 
formulated a realistic model for the LEC growth of axisymmetric 
crystals and have derived the differential equation relating the time 
evolution of the DWGS to radius and length. We show that the 
magnitude of the DWGS at the crystal's "shoulder" is inversely 
related to the radius of curvature. Furthermore, the meniscus by itself 
gives rise to a precursor or early warning in the signal, which means 
that the maximum in DWGS precedes the maximum in shape by a 
few hundred seconds. The existence of a secondary maximum or 
aftershock in the signal that is the sole consequence of liquid encap- 
sulation is also demonstrated. Excellent agreement has been obtained 
between DWGS and the signal predicted from the measured shape of 
a grown crystal. Thus, prospects for automatic diameter control are 
encouraging. 

I. INTRODUCTION 

GaAs is one of the key semiconductor materials that serves as a 
substrate for light-emitting diodes (LEDs), lasers, and field-effect 
transistors (FETs). Paralleling the development pattern of growth 
techniques for other single crystals, the issues of primary interest have 



477 



evolved from questions of quality (elimination of twinning and defects, 
doping uniformity) to that of economy of size. However, scaling up the 
dimensions of GaAs crystals grown by the liquid-encapsulated Czo- 
chralski (LEC) technique necessitates introducing sophisticated 
schemes for diameter control. By achieving satisfactory control, it will 
also be possible to run for extended periods of time with minimal 
supervision, and to attain higher crystal yields, as well as a reduction 
in defect generation brought on by shape change. 

Since diameter control in Czochralski growth has been the subject 
of an excellent recent review by Hurle, 1 here a brief outline of previous 
efforts with respect to LEC growth will suffice. Unlike Si 2 and a wide 
variety of oxide crystals 3 (e.g., GGG, LiTaOs) for which successful 
diameter measurement and control systems have been developed, the 
realization of a viable system for the LEC growth of III-V compounds 
has been much more difficult to achieve. This is largely due to effects 
associated with the growth chamber under high pressure, the presence 
of the encapsulating layer of B 2 O 3 (<0, and the phenomenon of anoma- 
lous density (d Uqu id > ckoiid) in some semiconductor materials. 

A number of approaches to the control problem have been taken. 
Pruett and Lien, 4 and van Dijk et al. 5 have employed an X-ray beam 
passing through the high-pressure growth apparatus for GaP, designed 
to make the melt a high absorber of radiation relative to all other 
materials in the radiative path of the system. The subsequent use of 
an image intensifier gave an accurate television picture of the growth 
interface. Small changes in diameter could then be electronically 
extracted from the video signal and utilized for its control. 5 

Alternatively, considerable success has been achieved meeting the 
demand for large, closely controlled diameter GaP single crystals using 
a floating die technique. Cole et al. 6 employed a Si3N 4 ring floating in 
the GaP melt beneath the B 2 O 3 («0. This modified form of Stepanov 
ring 6,7 creates a long-term stable growth regime, permitting diameter 
tolerances of ±1 mm for a 50-mm diameter boule. The technique has 
also been applied to GaAs with different degrees of success, apparently 
depending on the crystallographic growth direction. 8 

In a third technique, advocated by Bardsley et al., 9 a weight signal 
or a derivative weight-gain signal obtained by detecting the apparent 
weight change of either the crucible or crystal during Czochralski 
growth is used to control the power output and consequently the 
diameter. Bardsley et al. have considered both the theory 9 and its 
implementation by an analogue servo-system. 10 Although these au- 
thors have proposed reasonable initial postulates and expressions to 
calculate the weight gain, the detailed analysis is limited only to small 
perturbations in diameter and the formalism precludes growth by the 
LEC technique. 
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For some time now we have routinely measured the apparent weight 
of LEC-pulled GaAs crystals with a high sensitivity load cell placed in 
the high-pressure chamber in series with the crystal pull rod. The 
output of the cell is recorded as a dc level. Then, the derivative of the 
signal is taken electronically with a device developed especially for 
that purpose, and is also recorded. This activity has served as a useful 
qualitative guide in the manual control of diameter. The major objec- 
tive of this paper is to develop the fundamental theory that governs 
the relationship between the shape of an axisymmetric crystal with 
arbitrary variations in its diameter (including the shoulder) grown by 
the LEC technique and the instantaneously detected derivative 
weight-gain signal. Besides gaining new physical insight, we expect 
that by means of these investigations we can focus on the important 
physical factors and analytical techniques essential to the eventual 
control of crystal diameter. 

As a first step, a tractable model is formulated for LEC growth with 
a meniscus. The treatment leads to an explicit expression for the 
derivative weight-gain signal exclusively in terms of the crystal cross 
section and its first and second derivatives, in addition to geometrical 
and material parameters. Next, the numerical techniques required to 
perform computer simulations are outlined. Then, the key features of 
the signal are investigated by using the probability density function of 
the lognormal distribution for the crystal contour. Among the effects 
considered in detail are the shape of the shoulder and the presence or 
absence of the B2O3(<0 encapsulant and meniscus. Furthermore, the 
derivative weight-gain signal is predicted for four contour lines — gen- 
erated by consecutive 90° axial rotations — of a specially prepared 
GaAs boule and compared with the signal measured during growth. 
Finally, in light of these modeling calculations, the prospective tech- 
niques for diameter control are discussed. 



II. THEORY 

In this section, by a consideration of the relevant features of the 
LEC growth process, a realistic model is developed for the derivative 
weight-gain signal in a form suitable for computer analysis. To make 
the calculations tractable, but without a serious loss of generality, the 
following set of assumptions and constraints has been introduced: 

(i) The crystal and crucible are axisymmetric with circular cross 
sections. Hence, the crystal's radius, r, can always be prescribed as a 
function of the axial location, L, which is monotonically related to the 
time elapsed since seeding. 

(ii) The crucible containing the melt and the encapsulant is a 
right circular cylinder of radius, R. Less restrictively, departures from 
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a right cylinder are permitted below the liquid level prevailing at the 
completion of growth. 

(Hi) The effects of crucible and crystal rotation and the concomi- 
tant viscous drag are neglected. 

(iv) As the initial melt level falls by distance € (Fig. 1), the entire 
amount of the liquid is transferred to the solid. 

(v) Seepage from the B 2 3 (<fl encapsulating column of mass, jub, 
and volume, Vb, to the melt/crucible interface (walls) is minute. 

(vi) The solid-liquid interface is planar. This allows the separation 
of the associated heat transfer problem from the calculations. 

(vii) Capillary forces between the melt, B 2 3 (/), and solid result in 
a meniscus of height, h, with a subtended angle between the vertical 



INITIAL 
MELT LEVEL 



B 2 3 (l) 





— B 2 3 (l) 



Fig. 1— (a) Schematic geometry of LEC pulling, (b) Enlarged view of the solid-liquid 
interface region. 
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and the tangent to the "skirt" of the meniscus (joining angle), a (Fig. 

(viii) As the crystal is pulled at constant rate p, contact with the 
melt is uninterrupted. 

(ix) Physical constants required in the calculations are average 
values representative near the melting point of the crystal. 

The preceding conditions in combination with the geometrical de- 
scription of LEC pulling in Fig. 1 and Archimedes' principle permit 
writing the weight signal, W„ corresponding to detected crystal mass, 
m 8 , as the sum 



W. = m.g = 



W, 
true 

weight 



-Wi 

buoyancy 

correction 

to Wi 



+ W 3 

vertical 

projection 

of capillary 

force 



+ W< 

weight of the 

cylindrical 

meniscus 



-Wi 

buoyancy 

correction to 

W 3 + W< 



(1) 



where 



Wi = (m e + m g )g = mg 

W 2 = V eP Bg 

W3 = 2irra cos a 
W 4 = r 2 mhpig 

W 5 = V men pBg 



(2a) 
(2b) 
(2c) 
(2d) 
(2e) 



The symbols m, V, and p without subscripts designate the mass, 
volume, and density of the growing crystal, respectively. The same 
letters with subscripts B, % e, and g refer to the B 2 3 encapsulant, 
melt, and that portion of the crystal that resides either in the encap- 
sulant or the gas phase, respectively. The symbols a and Vm en represent 
the surface tension and meniscus volume, respectively. 

Since V e = (m — m g )/p, with the abbreviation k x = p B /p, W\ — W 2 
in eqs. (2a) and (2b) can be transformed into 



Wi - W 2 = m{\ - ki)g + kimgg. 



(3) 



Then, differentiating eqs. (3), (2c), (2d), and 2e with respect to time, t, 
and using eq. (1), we obtain the derivative weight-gain signal (DWGS) 
in the form 

dW a /g ,,dm , dm e 2tt dr 

— — — = (1 — k\) —— + k\ — — H a cos a — 

at dt dt g dt 

2mo . da ( ,dr , . dh\ dV men tMX 

__ rsma _ + p ^ 2 ^_ + r2 ,_J_ pB __. (4) 

It can be readily shown that dm/dt is a function of the other 
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variables. Clearly, if the crystal is axisymmetric, one has 



dm 2 dL 

-dF = p " r W 



(5) 



Moreover, the continuity of growth provides the relation 

L = f+pt-h, (6) 

or in derivative form 



dL = df _dh 
dt dt P dt' 



(7) 



In addition, the conservation of matter during crystallization as the 
liquid level falls leads to 



TO = pA.TR 2 t— Vmen) 



or 



/ D2 d«f dV mea \ 
\ wR dt dTJ- 



dm I D 2<M dVmen 

-dt =p '[ wR It dt 



(8) 



(9) 



A combination of eqs. (5), (7), and (9) permits the elimination of 
d€jdt\ hence eq. (5) can be rewritten as 



dm 

~dt 



= pnr 2 



(kdVmen 2 _ 2 ffi \ 

7T dt +Pr r ~dt) 



(10) 



where k = p*/p and R 2 = kR 2 - r 2 . 

Substituting eq. (10) into eq. (4) for the DWGS yields after several 
algebraic operations 

dW 8 /g dm g p(l - kfrrp£*R* 2tto da 

— — ^^— = «i —— t _„ r sin ft — — 

dt * dt R 2 g dt 



dV B 



dt 



R 2 

~2 



(1 ~ ki)p,-=i ~ PB 



2 dh 

+ **» 

+ 2it — 
dt 



R 2 



rnpi H — cos a: 
8 



. (ID 



In essence, the general form of the DWGS is given by eq. (11). Since 
all the variables directly depend on L or indirectly through r on L, it 
is convenient to replace the derivatives d/dt by d/dL-dL/dt. Accord- 
ingly, we have 
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dW 8 /g . dm g dL p(l - ki)irp,r 2 R 2 lira da dL 

-dT =ki -dL-dt + W '~ rmk °dL* 

dhdL 
R 



dVmen dL 



dL dt 



r 2 



(1 - kx)p,-=i - p B 



+ rw lLTt l " 



i-d-*i>5p 



dr dL ( , a \ tmM 

+ 2 "dL*\ rh >" + g COSa } <12) 

Equation (12) completely describes the DWGS once the few remaining 
unknown functions are evaluated. These are the axial derivative of the 
crystal mass outside the B 2 C>3(<0 (dm g /dL), the macroscopic growth 
rate (dL/dt), and the meniscus (h, dh/dL, V men , and dV men /dL) in 
terms of r(L) and the required physical and geometrical constants. 
The crystal shape r(L) we consider here as an input signal and the 
choice of different functional forms will be postponed to the next 
section. 

2. 1 Evaluation of dm„/dL 

From Fig. 1 we conclude that the total cylindrical volume occupied 
by the segment of the growing crystal in V e , the encapsulant, and the 
meniscus are conserved. Thus we can write the integral relation 

J r 2 irdL - | * r\dL + Vn™, + V B = wR 2 (L -L g + h), (13) 
Jo Jo 

where the two integrals on the left-hand side are equal to 

rL 

r 2 (L)7rdL. 



I 



Thus we can express m g /p as 

m g /p = J r 2 irdL = | r 2 irdL + V mea +V B - itR 2 (L -L g + h). 

Jo Jo 

Differentiating the above equation with respect to L gives 

dm g /p = r2 „ + dV^ L _ 7TR2 ( 1 + dh_ dLA (14) 



dL dL \ dL dL, 

To obtain dL g /dL we define a function <j>(L, L g ) =0 by subtracting 
the right-hand side from both sides of eq. (13). Then, employing the 
well-known result for the differentiation of implicit functions we have 

dL g _ d<f>/dL 
dL d<f>/dL g ' 
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where 



and 



^«r%L)* + dymen 



dL dL 



--(■♦S) 



^=-rHL g )Tr+*R 2 . 

The actual value of L g can only be determined by a numerical proce- 
dure. From eq. (13) we get 

f L r 2 TTdL + — +V mea | r 2 dL 

~ J 1 -(L + A)=A___ L „ (1 6) 

where the left-hand side includes only known quantities. 

While the entire crystal is submerged in B 2 3 , m g = and 
dnig/dL = 0. In that case the counterpart of the volume balance in eq. 
(13) is of the form 

r 2 *dL +V B + V m «n = irR 2 (L + h + A), (17) 

'0 



f 



where A is the distance of the tip of the crystal from the B 2 3 (/) — gas 
interface. Rearranging eq. (17) yields 



1 



L 

r 2 -ndL + — + Vu 



J^ (L + h) = A. (18) 



irR* 

Comparing eqs. (16) and (18) we observe that the left-hand sides are 
identical. As the crystal protrudes through the encapsulant, A changes 
from positive [eq. (18)] to negative [eq. (16)]. This provides an impor- 
tant clue in the computation of L g . The numerical solution of eq. (16) 
for L g starts when A = 0. 

2.2 Determination of the macroscopic growth rate, dL/dt 
The macroscopic growth rate given in eq. (7) can be rewritten as 
dL _ dtdL _dhdL . . 

~dt~dl~di dldt +P ' 



Hence, by a rearrangement of eq. (19) dL/dt becomes 

dL p 

~dt~ 1 - df/dL + dh/dL' 
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(20) 



To eliminate the melt coordinate £ we use the relation for the conser- 
vation of melt [eq. (8)], which gives 



Jo 



pr*dL + p/Vm«k 
,_ m + ptVum _ Jo 

psrR 2 p^R 2 • (21) 

Differentiating eq. (21) with respect to L yields 

rf^ = jf_ dV men /dL 

dL kR 2 + ttR 2 ' (22) 

Finally, introducing eq. (22) into eq. (20) gives the macroscopic growth 
rate in the form 

dL = p 

1 W ~dir/ ) ^ 

2.3 Transformation of the crystal length to time coordinate 

In general, r is known as a function of L. On the other hand the 
DWGS is recorded as a function of time. Therefore, a transformation 
from L to / is necessary in the analysis. From the continuity condition 
[eq. (6)] t can be expressed as 

t _ l + h(L)-AL) 

P 

Substituting eq. (21) into eq. (24) yields 

/. 
r 2 dL 

t= \ L + h ~ :2L w 2 — ^^l ' p - (25) 

2.4 Meniscus height and related properties 

Mika and Uelhoff n have determined the meniscus shape and inter- 
face height, h, occurring during Czochralski growth by a numerical 
solution of the Euler-Laplace differential equation. Although, in prin- 
ciple, numerical results for the meniscus could be generated concur- 
rently with the computation of the DWGS, here we prefer to rely on 
faster closed form solutions of the capillary equation. Fortunately, as 
shown by Mika and Uelhoff, 11 the analytical approximation of Tsivin- 
skii 12 describes the exact values of h with excellent precision over a 
wide range of joining angles, a, and radii. We have verified that for 
a > -10° and r > 0.5 cm the error in using Tsivinskii's approximation 
is less than ~2 percent. Subsequent computer simulations have dem- 
onstrated that in practical situations, except when there is a very rapid 
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drop in radius, the stated limit on a is obeyed. The fact that for a short 
time after seeding r is less than 0.5 cm is of minor importance owing 
to the initial insensitivity of the DWGS to h. 

Therefore, Tsivinskii's equations 12 for the interface height and shape 
will provide the basis of our DWGS analysis. In our notation Tsivin- 
skii's equation for h is of the form 

h = A[(l - sin a + u 2 ) y2 - u], 



(26) 



where 



u = 



A cos a 
4r(L) 



and A is a constant. According to Egorov, Tsivinskii and Zatulovskii's 
modification for LEC growth of Tsivinskii's original work A is given 
by 13 



A= y/^ 



2a 



PB)g 



(27) 



The joining angle a is the sum of the growing angle, <f>, and contact 
(wetting) angle, /8 (Fig. 1), i.e., 



a = <f> + p. 



(28) 



In general, /? is a constant, while the growing angle is a function of the 
derivative dr/dL according to 9 



dr 
tan* = — , 



* = tan->-. 



(29) 



Thus, h [eq. (26)] can be expressed in terms of r and dr/dL at the 
solid-liquid interface. 

Consequently, the required first derivative of h, dh/dL, includes the 
second derivative d 2 r/dL 2 . This can be readily shown by the substi- 
tution of the combined eqs. (28) and (29) into eq. (26). Then, using 
standard trigonometric relations for angle-sums and inverse functions 
and differentiating the expression thus obtained with respect to L we 
find 



dh 
dL 



rt). 



uv — X 



(1 - sin a + u 2 ) 1/z 



— v 



(30) 



where 



\ = 



1 d 2 r/dV 

2 (dr 

l+ (z 



n dr . n 
cos p — ==■ sm p 
dL 
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V= ~4? 



rd 2 r/dL 2 f dr 
— —^-cos^sxn 



mp\ 



dr ( n dr . 
+ dL[ COS/ *-dL 8mf t 

Another derivative of interest in eq. (12) is da/dL. From eqs. (28) 
and (29) we have 

da d(tan l dr/dL + P) 1 d 2 r 



2 jr2- 



dL dL 1 + (dr/dLf dL 



(31) 



A property closely related to the interface height is the volume of 
the skirt-shaped meniscus, V men (Fig. 1), and its axial derivative dV men / 
dL. The quantity V men can only be evaluated by a numerical technique. 
Egorov et al. 13 provide the x, y coordinates of the meniscus in the form 
of the integral 

*- I -j==dy + r(L), (32) 



where 



(1 . ,1 cos O I , , 



Owing to axial symmetry summing up narrow segments of area x 2 tt 
provides the volume of the meniscus as 

V me n = 7r5*£ x 2 ( yi ), (33) 

y,-° 

where 8 is the thickness of a segment. 

Apart from performing numerical differentiation, no simple method 
exists for the evaluation of dV mea /dL. If the shape of the meniscus is 
that of a right cylinder, then Vmen and dV men /dL are given by 

V cyl = r 2 irh 

and 

dVcyi dr 2 dh 

_ =2r7rA _ + r V_. (34) 

Simulations indicate that the numerical derivative dV men /dL is closely 
approximated by 

tf Vmen ^ V m en dV C y\ 

dL Vcyi dL 
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III. RESULTS AND DISCUSSION 
3.1 Computing methods 

We have succeeded in reducing the theoretical determination of the 
DWGS to a differential equation in r, dr/dL and d 2 r/dL 2 . Therefore, 
if the shape of the crystal is given, the DWGS is calculable; conversely, 
the DWGS completely defines the radius of the crystal at any moment 
of interest. Here, we address ourselves to the first part of the problem, 
though the eventual objective is diameter control emphasizing the 

DWGS. 

The program to evaluate the DWGS was written in HP Basic. A 
key feature in the efficient organization of the computations is the use 
of subroutines for dm g /dL [eq. (14)], dL g /dL [eq. (15)], dL/dt [eq. 
(23)], h [eq. (26)], dh/dL [eq. (30)], da/dL [eq. (31)], V men [eq. (33)], 
dV m en/dL [eq. (35)], as well as for r(L) and its first and second 
derivatives. 

Numerical integration was necessary to obtain V men and the crystal 
volumes between O and L and O and L g . It was found advantageous 
to evaluate the crystal volumes in the main routine because during the 
simulated growth of a crystal incrementing L in steps of 0.1 cm, the 
accumulating volume elements could be retained in memory. 

As mentioned earlier, the quantity L g is only required after the seed 
emerges from the B 2 (W) encapsulant. Mathematically speaking, this 
occurs when the left-hand side of eq. (16) or (18) becomes negative. To 
determine L g , at each increment of L the left-hand side of eq. (16) is 
evaluated. Then, a linear search is performed, substituting trial values 
of L g into the right-hand side of eq. (16) in the range L' g - 0.2 to L' g 
+ 0.4 cm in steps of 0.01 cm, where L' g is the previous solution to eq. 
(16) at L - 0.1 cm. In case the 0.01-cm mesh around L' g was too coarse 
to yield a solution, an option with 0.001-cm increments was available. 

There are several mathematical tools available to describe the shape 
of the grown crystal, r(L). From the actual cross section a large number 
of r, L coordinates can be generated in tabular form. Then, one obvious 
choice is to employ an interpretation scheme based on point-by-point 
tabular values of r(L) in addition to numerical first and second deriv- 
atives. A more sophisticated alternative is to determine the Fourier 
coefficients of the crystal shape. However, we have learned by expe- 
rience that the most efficient and satisfactory scheme to convert the 
tabular data into functional form is by using cubic spline regression. 

A cubic spline function is a piecewise polynomial of degree three 
that joins adjacent polynomials in "knots." At the knots the functional 
values as well as the first and second derivatives are continuous. 
Ordinarily, a spline passes through all the points within the interval 
bounded by two knots. In contrast, in spline regression a cubic least- 
square fit for the same points is obtained. We have adapted, for the 
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present application, Wold's spline regression technique in terms of B- 
splines. 14 

Fortunately, the key features of the DWGS can be illuminated 
without invoking advanced data-fitting techniques. A judiciously cho- 
sen elementary transcendental function representing r{L), which has 
n continuous derivatives, will be found suitable to illustrate the effect 
of various factors governing the course of the DWGS. 

3.2 Major influences on the derivative weight-gain signal 
3.2.1 Crystal shape 

Typical LEC crystals, including GaAs, exhibit a pronounced 
"shoulder" as the radius rapidly expands from the seed. Boule-to-boule 
variation in the rate of approach to and radius of curvature of the 
shoulder is not uncommon. A mathematical description of the shoul- 
der, in general, and a prototype crystal, in particular, is possible by 
means of the probability density function of the lognormal distribu- 
tion. 15 Accordingly, we can write for a seed diameter of 0.6 cm 

Qn 2 L/L m ) 



r(L) = C 



exp — 



.sLJfrr 2s 



+ 0.3. (36) 



The median length, L m , and the constant s are chosen such that the 
maximum r is 2.8 cm at L = 2.2 cm. From these conditions we have 

L m = 2.2 exp(s 2 ) 

and a lengthy expression for C that is also only dependent on the 
standard deviation, s. Hence, the function r(L) can be constructed if 
s is given. 

In Figs. 2, 3, and 4 we present the crystal shape r(L) for s = 0.4, 0.8, 
and 1.2, respectively. It can be seen that with decreasing s the radius 
of curvature at the maximum decreases (i.e., \d 2 r/dL 2 \ increases). 
Furthermore, it can be shown that the rate of approach to the shoulder 
(dr/dL) rapidly rises as s diminishes. 

Substituting eq. (36) and its first and second derivatives into the 
DWGS [eq. (12)] and the subsidiary equations yields the signal and 
the crystal cross section as a function of time, both of which are shown 
in Figs. 2, 3, and 4. The required parametric values are listed in Table 
I. The results demonstrate that the DWGS is extremely sensitive to 
dr/dL and | d 2 r/dL 2 1. In fact, the steeper the rise in the shape, the 
sharper the peak in the DWGS and the larger its absolute magnitude. 

Two other characteristic features of LEC pulling that are of crucial 
importance can also be observed in Figs. 2, 3, and 4. First, the maximum 
in DWGS precedes in time the maximum in shape by a few hundred 
seconds. This time lag we shall designate as the "precursor." The 
second property is the additional maximum in DWGS at a time when 
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Fig. 2— The DWGS and crystal cross section as a function of time for LEC pulling. 

The shape versus length curve ( ) is based on eq. (36), setting s = 0.4. The time 

derivative of the unencapsulated weight [dm g /dt, eqs. (14) and (19)] is also shown. 
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Fig. 3 — The DWGS and crystal cross section as a function of time for LEC pulling. 

The shape versus length curve ( ) is based on eq. (36), setting s = 0.8. The pull rate 

term in eq. (12) — the factor usually employed in conventional diameter control — is also 
shown. 

the crystal's radius has already declined. We shall refer to this phe- 
nomenon as the "aftershock." 

We have investigated the source of the precursor and aftershock. 
Simulations show that neither r nor its derivatives bears any respon- 
sibility. Examining the variation of the interface height with time (Fig. 
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Fig. 4 — The DWGS and crystal cross section as a function of time for LEC pulling. 

The shape versus length ( ) curve is based on eq. (36), setting 8 = 1.2. The meniscus 

height [eq. (26)] is also shown. 



Table I — Parameters for the calculation of the 
DWGS 

Density of crystalline GaAs, d[g/cm 3 ] 5.17 

Density of molten GaAs, d,[g/cm 3 ] 5.71 

Density of B 2 3 (^), d fl [g/cm*] 1.55 

Mass of B 2 O 3 («0. m B [g] 129.5 
Pull rate, p[cm/s] 3.6 X 10" 4 

Crucible radius, .R[cm] 3.9 

Wetting angle, /?[deg] 15 
Capillary constant, A[cm] 0.4 

Surface tension, o[dyne/cm] 333 



4) we note that starting with a small but finite value at the seed, h 
saturates beyond the shoulder (h = 0.38 cm).* Hence, no clue can be 
extracted from the form of h. Concentrating on the term including the 
pull rate in eq. (12), one concludes that this conventional description 
of the DWGS 3 holds quite well in the early phase of growth but departs 
from reality near the shoulder and beyond (Fig. 3). The maximum in 
the p term and the shoulder perfectly coincide and no secondary 
maximum appears. 

If, however, one examines the derivative dm g /dt [eqs. (14) and (23), 



* The fact that at t =- 0, r ¥> 0, thus h > is finite contradicts eq. (24). Therefore, to 
correct eq. (25) the residual time 



is always subtracted from p. 



h(L = 0) - 



V men (L = 0) 



wR 1 
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Fig. (2)], it is found that this quantity rises from zero to a maximum 
value precisely at the time of the aftershock and that the magnitude 
of the DWGS in this regime is essentially dm g /dt. In view of the fact 
that dm-g/dt is specifically associated with B 2 O 3 (<0 encapsulation, a 
reasonable hypothesis is that the aftershock is a consequence of the 
LEC technique. 

3.2.2 Meniscus shape and liquid encapsulation 

To isolate the factors leading to the precursor and aftershock we 
have evaluated the DWGS corresponding to the lognormal r(L) [eq. 
(36), s = 0.8] in the absence of liquid encapsulation. In this case the 
major equations still hold provided p B , Vj», ttib, m g , and dm g /dt are 
taken as zero. In Fig. 5 we present the DWGS with and without 
interposing a meniscus between the melt and the growing crystal. It is 
immediately obvious that in both instances the aftershock disappears 
from the DWGS. Moreover, if the meniscus is also removed, the 
precursor is missing, i.e., the times to reach the peak in DWGS and 
cross section exactly coincide. In fact by simultaneously eliminating 
B 2 3 (tf) and the meniscus from the growth system, only the contribu- 
tion of the conventional pull rate term 3 is retained in eq. (12). 
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Fig. 5— The DWGS and crystal cross section as a function of time in normal Czo- 
chralski pulling without encapsulation. Curves with and without a meniscus present are 

given. The shape versus length curve ( ) is based on eq. (36), setting s = 0.8. When 

the meniscus is absent, the pull rate term in eq. (12) exactly coincides with the curve 
shown. 
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Further confirmation of these effects is offered in Fig. 6. Here the 
calculation of the DWGS is repeated for LEC pulling without a 
meniscus. As we expected, though the precursor is missing the after- 
shock reappears. Apparently, up to and over the shoulder the DWGS 
is governed by thep term. Near the secondary maximum the derivative 
drrig/dt becomes the dominant factor. 

A more physical interpretation of the aftershock is possible by 
plotting the time-dependent radius, r(L g ), of the crystal as it just 
protrudes through the B2Oa(«0. In Fig. 6 we show both r(L) and r(L g ) 
as a function of time. Besides an expected displacement of the radius 
along the time axis one notes that the maximum in r(L g ) occurs exactly 
at the time of the secondary maximum in DWGS. From these consid- 
erations a surprisingly simple explanation emerges. As more of the 
shoulder region becomes uncovered from the encapsulant, the DWGS 
registers the rapidly increasing true weight gain as opposed to the 
previously measured apparent (buoyancy-reduced) weight gain of the 
crystal. Hence, at the time of the maximum aftershock both the freshly 
grown layers at L, submerged in B2O3(«0, and the sizable exposed 
portions at L g are detected. 

The discovery of the precursor and the aftershock has a significant 
bearing on the diameter control of GaAs. The precursor is an early 
warning signal that predicts a maximum in radius a few hundred 
seconds before the actual event. In other words the DWGs is already 
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Fig. 6 — The DWGS and crystal cross section as a function of time in LEC pulling, 
free of a meniscus. The radius is shown as a function of time at the solid-liquid interface, 
L, and at the top of the B 2 3 layer, L t . The pull rate term is also shown. The shape 
versus length curve ( ) is based on eq. (36), setting s =■= 0.8. 
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dropping while the diameter is still increasing. Thus the controller has 
sufficient time to react to an unwanted change. 

When the aftershock is observed there is a good chance that the 
controller misinterprets it as an unexpected gain in diameter and, if 
undesirable, will respond accordingly. Then, of course, a visible reduc- 
tion in diameter will occur. By the time this is noticed, corrective 
action taken cannot restore the loss in diameter control. At best, 
beyond the shoulder a "sinusoidal" cross section with small amplitude 
results. 

In view of these observations one cannot base the diameter control 
algorithm for GaAs on the customary p term 3 in eq. (12). To be sure, 
from the seed up to near the shoulder perhaps it suffices. However, for 
the bulk of the boule, control by means of such a traditional method 
is clearly inadequate. 

3.3 Comparison with experiment 

We have extracted the salient features of the DWGS arising during 
LEC pulling by the use of the lognormal probability function for the 
crystal's cross section. To compare the experimentally determined 
DWGS with the theoretical one the idealized crystal shape must be 
replaced by that of a GaAs boule grown by the LEC method under 
closely controlled conditions. Some of the growth parameters of inter- 
est for this specially prepared crystal are listed in Table I. The weight 
and length of the crystal are 523.3 grams and 6.6 cm, respectively. 

The crystal coordinates (r versus L) can be derived from a two- 
dimensional projection employing a digitizer. Owing to the eccentricity 
of the actual < 100) boule, a single view is insufficient. We have obtained 
planar projections with sharp contours by a photographic technique. 
The crystal was placed on the top of a light box, backlit, and photo- 
graphed using a high-contrast film. Then, the procedure was repeated 
following a 90-degree rotation around the axis. In this manner, the two 
photos provided a total of four cross sections. Digitizing the shape in 
0.5-mm intervals 133 pairs of (r, L) coordinates were collected for each 
of the four contours. At the shoulder the maximum error between the 
crystal and its projection is less than 0.5 percent. The computed mean 
weight is 528.6 g with a standard deviation of ±9.3 g, which should be 
compared with the measured weight of 523.3 g. Clearly, an excellent 
tabular description of the crystal's cross section has been accomplished. 
The data points for one of the four views are shown in Fig. 7. 

The alternatives to represent the tabular data have been outlined 
earlier. In principle, because of close spacing, there would be no 
obstacle to a simple interpolation of the r(L) values. However, as we 
illustrated in Fig. 7, using the numerical first and second derivatives 
may be problematic on account of noise. Indeed, when the DWGS was 
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Fig. 7 — One of the four experimental contours of a Ga As crystal and its cubic spline- 
regression representation (— ). The cubic spline regression of the numerical first deriv- 
ative dr/dL is also shown. In addition, the numerical second derivative d 2 r/dL 2 is given, 
together with the quadratic expression derived from the cubic spline of dr/dL. 



computed from the numerical derivatives severe discontinuities ap- 
peared in the results. This suggests that the steep "jumps" in d 2 r/dL 2 
seriously affect at some locations the shape and magnitude of the 
DWGS. 

A different approach is the finite Fourier transform of the r(L) 
points. However, owing to slow convergence a large number of terms 
is required. Eight terms were found to be sufficient to describe the 
seed and shoulder regions. Beyond the shoulder 18 terms were neces- 
sary to fit the data, albeit with concomitant high-frequency noise. 
Since the number of required terms was unpredictable and the com- 
putation inefficient, Fourier representation of the r(L) data was aban- 
doned. Nonetheless, the DWGS based on Fourier analysis provided a 
reasonable description of the measurements. 

Cubic spline regression 14 has proven to be the most promising 
method to cast the tabular data on crystal cross sections into functional 
form. At the middle of the axis, the crystal was separated into two 
slightly overlapping domains and within each 12 equally spaced knots 
were positioned. The outstanding spline fit to one of the contours is 
given in Fig. 7. There are two avenues open to obtain the first 
derivative dr/dL. One is a simple differentiation of the shape's spline, 
which results in a continuous quadratic representation. Unfortunately, 
in this case the second derivative becomes a continuous network of 
connected straight lines (the third derivative is discontinuous) leading 
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to a jagged DWGS. A more fruitful approach is to fit the numerical 
first derivative itself by a cubic spline. Then, the noisy second deriva- 
tive is described by a quadratic expression. In Fig. 7 we show the cubic 
spline fit to the numerical first derivative data as well as to the 
quadratic function passing through the numerical d 2 r/dL 2 . We note 
that a reasonably good description of the derivatives has been obtained 
by the spline regression technique. Therefore, all four contours were 
treated in a like manner. 

Having thus established the shape and shape derivatives for an 
actual GaAs crystal, one can readily evaluate the corresponding DWGS 
by means of the analytical tools described earlier. Among the input 
parameters listed in Table I, all the densities were obtained from a 
recent critical evaluation of Jordan. 16 The selected contact angle (15 
degrees) is consistent with theoretical and experimental findings on 
Ge and Si. 17,1819 The capillary constant A and the surface tension, a, 
are connected via eq. (27). Though the recommended value for A by 
Egorov et al. is 0.48 cm, 13 we have achieved a slightly better fit to the 
experimental DWGS using A = 0.4 cm, which is equivalent to a = 333 
dynes/cm. The surface tension of GaAs thus obtained is consistent 
with that for other III-V compounds (InP 20 and GaSb 21 )Si and Ge. 22 

The gross appearance of the DWGS is not overly sensitive to a 
change in A. Nonetheless, a detailed examination reveals that the peak 
corresponding to the shoulder becomes higher and broader as A 
increases from 0.28 to 0.48 cm. With respect to the rest of the DWGS 
one can conclude that the peak-to-valley dynamic range diminishes as 
A drops. Among the effects of the other parameters, we have observed 
a rise in DWGS with an increase inp and a reduction of the time scale 
with a decrease in crucible radius. 

In Fig. 8 we have reproduced the experimental derivative weight- 
gain signal as a function of time, as obtained on a carefully calibrated 
strip-chart recorder during the LEC growth of GaAs. The predicted 
DWGSs for the four contour lines as well as the time-dependent cross 
sections are also shown. Clearly the theoretical curves envelop the 
measured values, providing an excellent overall description. 

The computed DWGSs exhibit both the early warning precursor 
and aftershock. In the inflection before the first peak we can discern 
the time when the seed first protrudes through the B 2 3 (/). Since an 
actual crystal may possess additional cross-sectional bulges beside the 
shoulder, the secondary maximum in DWGS reflects both the memory 
of that earlier maximum in radius and the current expansion. To a 
small extent the DWGS is influenced by the input shape derivative 
function. For example, in the trough above the shoulder, a spline of 
the first derivative with many more knots would lead to a smoother 
transition to the next peak. Likewise, near the last peak, an improved 
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Fig. 8 — DWGS ( — ) and crystal cross section as a function of time for four experi- 
mental contour lines of LEC-pulled GaAs. The measured DWGS was traced from the 
original strip-chart recording. 

agreement with the experimental data is possible at the cost of a finer 
derivative spline. Near the bottom of the crystal the accord with the 
measured DWGS is more limited, perhaps on account of the increas- 
ingly negative turn in the joining angle a (between -20 and -40 
degrees). Under those circumstances Tsivinskii's approximation for 
the meniscus properties becomes inaccurate. 11 

3.4 Prospects for diameter control 

The excellent agreement between predicted and experimental 
DWGS strongly suggests that diameter control of GaAs by the deriv- 
ative measurement is feasible. There are two approaches to diameter 
control implied by the modeling calculations. One is a real-time in- 
stantaneous determination of the radius during growth from the 
DWGS by solving the differential equation (12). Then, depending on 
whether r is increasing or decreasing from a preset limit, the RF input 
power is changed by a suitable amount. 

The other method is the construction of an "ideal" crystal with 
appropriate tolerances on a drawing board. Taking into account the 
tolerances for this crystal, a band of DWGS can be evaluated as has 
been done for Fig. 8. Subsequently, the control of diameter is reduced 
to the problem of adjusting the power input in such a fashion that a 
template established by the band of DWGSs is followed. 
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